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Abstract 

I describe how nuclear structure can be predicted from lattice QCD through low- 
energy effective field theories, using as an example a world simulation with relatively 
heavy up and down quarks. 


1 Introduction 


Weakly bound systems are fascinating for the surprising quantum features they display 
regardless of the details of their short-distance structure. They beg for a description with 
effective held theory (EFT), because that is the general framework to turn separated 
physics scales into a controlled expansion based on symmetries, rather than details of 
the dynamics. I have been smitten with EFT since shortly after leaving Sao Paulo for 
my Ph.D. in Austin. Only much later did I learn that a master of weakly bound systems 
(halo nuclei |1], neutron-rich nuclei [2], atoms near Feshbach resonances [3],...) lived in my 
hometown. This contribution is dedicated to this master. Professor Mahir Saleh Hussein, 
on the occasion of his 70th birthday. 

I will argue here, based on the work of Refs. [HEIIG], that recent lattice QCD (LQCD) 
data m El El do] suggest that light nuclei are weakly bound even in a world with relatively 
large quark masses. It seems that these data can be described with an EFT, Pionless or 
Contact EFT [TTl fT8] , where all degrees of freedom except nucleons are implicit [H] . 
This is, mutatis mutandis, the same theory that one can use to describe other systems 
and processes characterized by relatively large sizes: halo nuclei such as ®He [15], shallow 
molecules such as He 4 trimers [T6|, atom recombination near Feshbach resonances ca. 
etc. It has been known for some time to apply to light nuclei at the physical value of the 
quark masses [HI Il9l EQl EH EH ES] . If more or less familiar nuclear structure can be 
predicted for larger quark masses, perhaps nuclear physics is less accidental than we are 
used to think. 

This is not just an academic exercise. As quark masses decrease in LQCD, pions 
become lighter and another, less universal but more predictive EFT, Pionful or Chiral 
EFT [m [2l| , can be used to connect results at different quark masses [25] . We will at 
that point be able to predict even more deeply bound nuclei from LQCD, following the 
same steps as in Refs. [lElE], just with Pionless EFT replaced by Chiral EFT. 

After describing the lattice world in Sec. [2] and reviewing Pionless EFT in Sec. [3l the 
main results for lattice nuclei are summarized in Sec. HI An outlook is offered in Sec. [5l 

1 hope Hussein is pleased with another unexpected gift offered by weakly bound systems. 

2 EFT and QCD 

It is intuitively clear, and supported by our experience in physics, that only certain de¬ 
grees of freedom and symmetries are relevant at a given distance scale. Less obvious, but 
equally supported by the existence of virtual processes in quantum mechanics, is that 
all interactions allowed by symmetry take place among the relevant degrees of freedom. 
EFT is simply the framework that incorporates these facts. To be predictive, the basic 
assumption (“naturalness”) is that once a few scales are identihed, the inhnite number of 
interaction strengths, or “low-energy constants” (LECs), can be written as combinations 
of these scales, times numbers of 0{1). As a consequence, observables, which can be 
expressed through T matrices for various processes, can be obtained as controlled expan¬ 
sions in Q/M, where Q represents the momenta of interest and comparable mass scales. 
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and M, high mass scales. (I use natural units where h = c = 1.) Because all interactions 
are included, observables are renormalization-group (RG) invariant, that is, independent 
of the arbitrary regularization procedure used to separate explicit from implicit degrees 
of freedom. For a basic introduction to EFT, see Ref. [26] . 

Nuclear physics is well described in terms of nucleons subject to (possibly approxi¬ 
mate) Lorentz, (possibly approximate) baryon-number, (approximate) time-reversal, and 
(approximate) parity invariance. QCD for the two lightest quark flavors, as relevant for 
nuclear physics, has essentially two separate scales: 

1- Mqcd 1 GeV, where the coupling constant of QGD formulated in terms of quarks 
and gluons becomes large. Mqcd sets the scale for hadronic masses, including the 
nucleon mass m^v, and for 4:71 where f.^ is the radius of the “chiral circle” formed 
by the set of SU{2)l x SU{2)ji ~ ^(^(d) minima of the QGD effective potential. 
Picking one of these minima leads to spontaneous symmetry breaking, the emergence 
of pions as Goldstone bosons, and the manifestation of ~ 90 MeV as the pion 
decay constant. 

2. m ~ 5 MeV, the average quark mass, which explicitly breaks chiral symmetry, 
creates an absolute minimum of the QGD effective potential, and endows pions 
with a mass = 0 {Mqcdt^)- Ft also affects most other quantities, including rriN- 
We can trade ffi for 

Gontrary to real experiments, LQGD simulations can probe worlds where takes 
different values. In fact, the high cost of light quarks roaming the lattice constrains 
present calculations to large values of While it certainly is a disadvantage that 
LQGD cannot reach realistic values yet, one can turn a disadvantage into an advantage 
by learning how nuclear physics depends on 

Existing laboratory and LQGD data for nucleon masses and light-nuclear binding ener¬ 
gies are summarized in Tab. [T] I also list the EFT results described in Secs. |3]and|ll The 
LQGD calculations are performed with equal light quark masses and no photons, hence 
isospin symmetry is exact and the binding energy Ba{S, /) is determined by the nucleon 
number A and the spin-isospin combination S, I of the state. It did not have to be so, 
but masses and binding energies in Tab. [1] increase more or less monotonically with m^r. 
Note that Ref. |2^ hnds no bound states in a large range of pion masses that includes the 
values in Tab. [H Since the model-independence of some of its results, obtained through 
a non-observable potential, remains a question mark, I do not consider these data here. 

From the experimental and LQGD data in Tab. (H we can infer the relevant momentum 
scales for real and lattice nuclei, which I list in Tab. |2l Besides pion and nucleon masses, 
the momentum associated with the excitation of the lowest baryon, the Delta isobar, in 
two-nucleon scattering [28] is also given, using the LQGD data for the Delta mass 
compiled in Ref. [2H]. As one can see, for the three values of the pion mass the typical 
nuclear momenta ^j2mNBA/A for A = 2, 3,4 are much smaller than m^v, suggesting that 
a description in terms of non-relativistic nucleons is always appropriate. 

EFTs for non-relativistic nucleons have some simple, useful features. Pair creation is a 
short-range effect and the theory can be formulated in terms of Pauli spinors representing 
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Table 1: Neutron and proton masses, and binding energies of the lightest nuclei at various 
values of the pion mass. All entries are in MeV. The hrst column summarizes experimental 
data, the third [10], fourth [7] and hfth [8] columns give LQCD data, and the second [23] 
and sixth [I] columns show EFT input (marked with *) and results. The EFT calculations 
are discussed in Secs. El and 01 
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Table 2: Momentum scales for various quark masses: nucleon mass, effective momentum 
for Delta isobar excitation, pion mass, and effective binding momenta for S'-shell nuclei. 
All entries are in MeV. 
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forward propagation in time. The EFT Lagrangian contains the usual non-relativistic 
kinetic terms in lowest order, with relativistic corrections implemented at higher orders 
in a Q/mj^ expansion. Nucleon energies are of 0{Q‘^/rriN)- Loops with antinucleons 
never need to be considered explicitly, and an V-body force does not affect the n-body 
system for n < N. At the V-nucleon level, all possible operators involving up to 2N 
nucleon helds are included, with an increasing number of derivatives. These interactions 
are highly singular and require regularization via some ultraviolet (UV) cutoff A. The 
explicit dependence on non-negative powers of A coming from loops is eliminated by 
renormalization of the LECs. Once (and only once) this is done at a given order, an 
integration over momenta in intermediate states contributes a factor 0 {Q^/4:7 t) to the T 
matrix. These factors of Q, rriisi and dvr, together with the sizes of the renormalized LECs, 
are the ingredients to build the Q/M expansion. 
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In nature, there is a large separation of scales, m.„ <C Mqcd- A nucleon thus consists 
of an outer cloud of pions at distances ~ l/m.^ surrounding an unresolved, dense core 
of size ~ 1/Mqcd- Empirically, nuclear sizes scale as Ra ~ where tq ~ 1.2 

fm is possibly set by a combination of and Nuclei are expected to be large 

on the 1/Mqcd distance scale because of the two effects: the size of the pion cloud 
around each nucleon and the piling up of nucleons. LQCD has to hght both effects to 
contain nuclei within the lattice length L > Ra, while striving for a much smaller short- 
distance regulator in the form of a lattice spacing 6 < 1/Mqcd- EFT offers a strategy 
to extrapolate QCD to the large distances involved in nuclear physics: i) calculate with 
LQCD y4-nucleon observables for A = 2,3, 4; ii) calculate the same observables with EFT 
and match LQCD, thus determining the LECs; and Hi) solve the EFT for A> 5 using the 
powerful initio” methods that have been developed in recent years, such as the no¬ 
core shell model (NCSM) |30], the effective-interaction hyperspherical harmonics (EIHIl) 
[31], and the auxiliary-held diffusion Monte Carlo (AFDMC) [32] methods. 

For momenta Q ~ m.,^ -C Mqcd, one can formulate an EFT, Chiral EFT [IHISI], which 
includes, in addition to nucleons, also pions and the lowest nucleon excitations. In this 
EFT, M ~ Mqcd, and one treats the inner nucleon cloud in a multipole-type expansion. 
The approximate chiral symmetry of QCD plays a crucial role, because it ensures that 
pions couple weakly at low momenta, which gives rise to a loop, or equivalently density, 
expansion for the larger, more sparse pion cloud. Chiral EFT allows one to, in principle, 
calculate the dependence of low-energy nuclear observables on [25]. Unfortunately, 
however, an RG-invariant formulation of Chiral EFT is still work in progress [33] . This is 
not too serious a problem in the sense that it is doubtful that the Chiral EFT expansion 
holds at pion masses explored so far by LQCD. For example, studies of the convergence 
of Chiral EFT for A = 0 suggest a breakdown of the expansion at a pion mass no larger 
than 500 MeV [3l]. When smaller pion masses can be reached and a proper formulation 
well developed. Chiral EFT could be used as a tool for extrapolation of nuclear quantities 
in as it is already for meson and one-nucleon observables. Chiral EFT LECs would 
then be determined from LQCD instead of experimental data, and a solution of Chiral 
EFT used to extrapolate LQCD to larger A. 

3 Pionless EFT 

In fact. Tab. |2] suggests that the extrapolation to larger A can already be performed. For 
all available pion masses, the typical nuclear momentum is not only much smaller than 
m^v, it is also smaller than m.^.. This is already true in nature and, as the pion mass 
increases, pion effects become more short-ranged relative to nuclear distances. Assuming 
that the Delta continues to be the lowest baryon excitation, any effects from other hadrons 
are also short-ranged. 

Thus, at momenta Q -C nucleons should suffice as explicit degrees of freedom. In 
the appropriate EFT, Pionless EFT |T3], M ~ and all interactions are of contact type, 
which in coordinate space give (renormalized versions of) delta functions and derivatives. 
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In nature, Pionless EFT describes well the properties of low-energy scattering and bound 
states for A = 2 [IB], 3 [H] 120], 4 jBH |22], and even (but less well) 6 [20] • One expects 
Pionless EFT to breakdown at some point as nuclei get denser, but its reach is presently 
unknown. 

In the two-nucleon sector, there are two independent non-derivative contact interac¬ 
tions, with renormalized LECs Ooi = C>( 47 ra 2 ( 0 , l)/m 7 v) and Oio = C>( 47 ra 2 (l, 0 )/m 7 v), in 
terms of the ^S'o and scattering lengths, respectively 02 ( 0 ,!) and 02 ( 1 , 0 ). Together 
with the general estimates for nucleon energies and loop integrals given above, each it¬ 
eration of this potential in the T matrix yields a factor 0{Q\a2\)- When |a 2 | 3> 1/M as 
a consequence of a shallow pole in the T matrix, one needs to include all iterations for 
Q> 1/1 02 1, namely solve the corresponding Schrodinger equation exactly. It is easy to 
show that the dangerous UV regulator dependence can be eliminated if the bare LECs 
C 5 /(A) oc 1/A. In coordinate space one can understand this by noticing that a delta 
function is a A^/dvr singularity and overwhelms the kinetic term, which grows at most as 
A^/miv, unless the associated LEG goes as 47 r/(mArA). The relative (o 2 A)“^ corrections 
in 6 * 57 (A) provide then just the balance necessary for a low-energy (real or virtual) bound 
state with binding momentum 0{l/a2)- The regularization procedure leaves behind a rel¬ 
ative error of 0{Q/h), which can be made arbitrarily small by taking A arbitrarily large. 
In leading order (LO) we obtain the hrst term in the effective-range expansion (ERE) of 
the T matrix. The residual cutoff dependence can be removed by the two-derivative inter¬ 
actions present in the same channels, which, as a consequence, have a natural relative size 
0{Q/M) with respect to LO. That means they make up the next-to-leading order (NLO), 
their renormalized LECs being of 0{4:Tia\r2/mi<4), with effective ranges |r 2 | = 0{1/M). 
The argument can be generalized to more-derivative terms in the S waves, which con¬ 
tribute at progressively higher orders. In higher waves, where there seem to be no shallow 
poles, all LECs scale with M according to their canonical dimension, which means that no 
other wave needs to be considered up to N^LO. RG invariance requires that, as allowed 
by their small relative sizes, subleading orders be treated in perturbation theory. The 
gory details are spelled out in Ref. m- 

In this way, Pionless EFT generates an expansion of the two-nucleon amplitude equiv¬ 
alent m to the ERE. At LO, only the two non-derivative contact interactions need to be 
included. Each LEG is determined by one datum, say the scattering length or the binding 
momentum of the shallow pole (which differ only by higher-order terms). At NLO the 
two S-wave two-derivative contact interactions need to be included in hrst distorted-wave 
Born approximation, and each new LEG requires another datum, say the effective range. 
A nice description of two-nucleon scattering at physical quark mass was found in Ref. 
[TB] . In the regime l/|a 2 | "C Q -C M, the amplitude is approximately scale invariant and 
SU{4:) spin-isospin symmetric [SB] . 

Nuclear-physics folklore would suggest that few-nucleon forces are of higher order. 
Indeed, in certain channels such as ^S' 3/2 neutron-deuteron scattering, high accuracy can 
be obtained in the hrst orders [2Q]. However, in the ^S' 1/2 channel [19] (and for three bosons 
in relative S waves [12]), the only way to eliminate non-negative powers of A (in particular, 
the “Thomas collapse” [3S] of the ground state) is to ensure that the non-derivative six- 
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nucleon operator (there is only one) is present in LO. Again the cutoff dependence of 
its bare LEG, -D(A) oc 1/A^, can be understood by a simple coordinate-space argument. 
When going from two to three bodies, the delicate balance between kinetic terms and two- 
body contact is destroyed because the number of kinetic terms doubles while the number 
of pair-wise interactions triples. The result is the Thomas collapse where the two-body 
attraction wins and leads to a three-body binding energy that grows as A'^/mN. This 
growth has to be canceled by the contact three-body interaction; since it involves two 
delta functions, or A®/(47r)^, the bare LEG should be roughly (47r)^/(m7vA'^). The cutoff 
dependence of D is still more complicated, though. At very low cutoffs the three-body 
force might be attractive or repulsive, but once A^/m^v exceeds the binding energy of the 
ground state, the three-body force must be increasingly repulsive to prevent the collapse. 
With a regulator procedure that preserves the approximate scale invariance of the two- 
body subsystems, at a critical A one can maintain the binding energy hxed only at the cost 
of making the three-body force attractive and accreting a very deep bound state. As the 
cutoff increases past critical, the two-body attraction continues to increase and the three- 
body force gets less attractive, till it becomes repulsive and a new cycle begins. In fact 
da EH] the bare LEG is on an RG limit cycle, approximate scale invariance is reduced to 
an approximate discrete scale invariance, and there is a tower of approximately geometric 
three-body “Ehmov states” [37]. Varying the renormalized D = Cl((47r)^a2/'mjv) shifts the 
position of the tower and changes the three-body scattering length, leading to a correlation 
known as the Phillips line [38]. In nuclear physics the explicit breaking of scale invariance 
is such that the “tower” consists of a single state, the triton. The two-derivative three- 
body force hrst appears at N^LO [IS] |19]. 

Galculations da [21], which are however somewhat limited in cutoff variation, indicate 
that four-body observables do not display non-negative powers of A up to NLO, in the 
absence of four-body forces. This is perhaps not surprising since the three-body force is 
effectively repulsive and the number of triplets grows faster than doublets. With hxed two- 
nucleon input, variation in D leads to a correlation between four- and three-body binding 
energies, the “Tjon line” [39], which passes close to the experimental point. Pionless 
EFT thus successfully postdicts the alpha-particle binding energy [2l] , and scattering can 
be calculated as well [22]. Notice, however, that the arguments above leave open the 
possibility of cutoff dependence in the regions where the three-body force is attractive, 
and indeed Ref. found sensitivity in four-body properties to a four-body scale. As 
far as I can see, there is yet no compelling argument that this sensitivity is due to an 
LO (or even NLO) force. I will assume that the dominant four-body force, presumably 
coming from the single non-derivative eight-nucleon operator, hrst contributes beyond 
NLO. Since, thanks to the Pauli principle, hve- or more-body forces involve at least two 
derivatives, they are likely of even higher order. 

It is probably safe to assume that, up to NLO, Pionless EFT is renormalizable with, 
besides two-body forces, a single non-derivative three-body force. Its LEG is determined 
by one three-body datum, say the neutron-deuteron ^S' 1/2 scattering length or the triton 
binding energy. Thus at LO (NLO) three (hve) data are needed as input, in addition to 
the nucleon mass, and everything else is a prediction. Pionless EFT is not just the ERE; 
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it is the extension to few-body systems that preserves model independence. Pionless EFT 
acconnts for a series of apparently nnrelated, qnalitatively nniqne phenomena, snch as 
the Thomas collapse, Ehmov states, the Phillips and Tjon lines, and presnmably similar 
correlations for bigger systems. With a few adaptations, it applies to other systems 
characterized by a small ratio r 2 /a 2 —see, e.g. Refs. [ISl UHl H^. The approximate 
discrete scale invariance has striking conseqnences for the spectrnm of few-boson systems 
im, where a state in Ehmov’s three-body tower generates pairs of “image” states in bigger 
systems. The alpha-particle gronnd and excited states can likely be interpreted this way. 
Moreover, since the LO three-body force is SU{4) symmetric, Pionless EFT provides a 
justihcation [121 [19] for the approximate S't/(4) symmetry proposed by Wigner [42] . 

It is nnfortnnately still nnclear how far np the nuclear chart this EFT can be pushed. 
The only calculation [23] beyond the four-nucleon system was based on the NCSM IMI, 
when the EFT Hamiltonian was diagonalized in a harmonic-oscillator basis. This basis has 
a natural UV cutoff in the form of a maximum allowed number of shells. It contains also 
an infrared (IR) cutoff provided by the spacing between shells. In the simplest approach, 
the LECs are htted to the experimental binding energies of the lightest nuclei for every 
cutoff pair, and binding energies for larger nuclei are calculated and extrapolated to large 
UV and small IR cutoffs. Results of an LO calculation [23], where the deuteron, triton, 
and alpha-particle ground-state energies were used as input in addition to the nucleon 
mass, are shown in Tab. [H The input data are indicated by a in Tab. [1] Estimating 
the error as 30% from r 2 /a 2 in the ^Si channel, one sees that the excited state of the alpha 
particle is postdicted very well, while ®Li is barely consistent. However, the error could 
be as large as 80% if we consider the ratio in Tab. [2] between alpha-particle momentum 
and pion mass. Higher-order calculations are clearly needed. 


4 EFT for Lattice Nuclei 

As noted in Ref. [1], the widening gap shown in Tab. |2] between pion mass and typical 
nuclear momentum implies Pionless EFT should work better at larger pion masses. In 
the hrst calculation ever to ht lattice nuclear data. Ref. [4] used the nucleon mass and 
light-nuclear binding energies at the highest pion-mass value, mjr = 805 MeV, from the 
NPLQCD collaboration [8] as input for Pionless EFT in LO. A calculation using the 
ruTT = 510 MeV binding energies from Ref. [7] is in progress [5]. The existence of a 
dineutron bound state allows the use of its binding energy instead of the alpha-particle’s 
as input. 

From Tab. [2] a conservative estimate for the Q/tUtt expansion parameter at = 805 
MeV is 40%. Two-nucleon scattering lengths and effective ranges are also available at 
this pion mass [9], and are consistent with an almost degenerate double bound-state pole 
in the T matrix of each S wave, which is thought to be incompatible with a short-range 
non-relativistic potential [43]. References 015] assume |r 2 | = 0{l/m.,^), as for physical 
quark masses. If one uses the ratio r 2 /a 2 , the error estimate for LO is instead 50%. A 
test of convergence will have to await an NLO calculation. 
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In a renormalizable theory, only convenience guides the choice of regulator. In the 
present case, we want, as in Ref. |11], a local potential that allows the use of many- 
body techniques that cannot handle non-local interactions well, such as AFDMC [32]. 
This can be achieved with a regulator function in the momentum transfer q, 

or its Fourier transform F(r^A^) in terms of the radial coordinate r. Ref. [4] employed 
two forms, fn{x) = exp(—with n = 1,2, which get increasingly closer to a sharp 
regulator. The isospin-symmetric Hamiltonian can be written in coordinate space as 

^ = -^5^V? + i^|3Gi„(A) + G„,(A) + (C,„(A)-C„,(A))<Ti.if,|nryA") 

^ i i<j 

+ + • • • ’ ( 1 ) 
i<j<k eye 


where (Tj/2 (Ti/2) is the spin (isospin) of nucleon i, J2cyc stands for the cyclic permutation 
of a particle triplet {ijk), and “... ” for terms containing more derivatives and/or more- 
body forces. The LECs rriN, C'io(A), C'oi(A), D{A), etc. depend on m^, since pions are 
part of the short-distance physics not included explicitly. The dependence of two- and 
three-nucleon observables in Pionless EFT has been studied with input from Chiral EFT 
in Ref. [45] . 

The two-nucleon Schrodinger equation was solved [1] for the LO Hamiltonian with the 
Numerov method, and C'io(A) and Coi(A) htted to the deuteron i?2(l,0) and dineutron 
R2(0,1) binding energies [8], respectively. The cutoff dependence of the LECs is found 
to be qualitatively similar to other regulators mm- Csi{A)A approaches a regulator- 
specihe constant at a rate determined by y/mJv^S/OS^T). For large cutoffs one should have 
in LO 02 ( 1 , 0) ~ 1/A/mAri? 2 (l, 0) ~ 1.12 fm. For cutoff variation in the range 2-14 fm“^. 
Ref. [1] hnds explicitly 02 ( 1 , 0 ) = (1.2 ± 0.5) fm ((1.1 ± 0.1) fm) with the regulator /i 
(/ 2 ). For comparison, LQCD gives 02 ( 1 , 0 ) = 1-821q;]31q;] 2 fm directly [9]. The situation 
is similar in the ^So channel. 

For systems with 3 < A < 6 nucleons the Schrodinger equation was solved [1] with 
the EIHH method, where the wavefunction is expanded into a set of antisymmetrized 
hyperspherical-harmonic spin-isospin states. Convergence is controlled by the hyper- 
angular quantum number Kmax, results being obtained by extrapolation to the limit 
Kmax —t C )0 [31]. The corresponding error was estimated to be smaller (for the lighter 
systems, much smaller) than the EFT truncation error. For systems with A > 4 the 
AFDMC method was also used. In this technique [52], the ground-state energies are 
projected from an arbitrary initial state by means of a stochastic imaginary-time prop¬ 
agation. The numerical simulations are simplihed by the introduction of auxiliary helds 
via a Hubbard-Stratonovich transformation. In all these calculations [1], the regulator 
employed was /i with 2 < A fm < 8. 

The LEC D{A) was determined [4] imposing that the ^H/^He binding energy is 
reproduced at any value of A. It was found that D(A)A^ approaches a hnite limit, as for 
other regulators [T21 [T9] . The limit-cycle behavior is not seen, as in other cases when the 
number of three-body bound states is kept hxed with regulators that do not preserve the 





approximate scale invariance of the two-body subsystems, e.g. Ref. [IB]. With LECs thus 
hxed, a complete LO potential is available to predict other properties of lattice nuclei. 

The four-nucleon system, solved [I] with both EIHH and AFDMC methods, provides 
a consistency check between the two ab initio methods, and between them and LQCD. 
The two ab initio methods produced results for the ^He binding energy R 4 that agree well 
within the (large) LQCD error. In either case, B 4 was found to depend only weakly on 
the cutoff, changing by about 20% when A grows by a factor of 4. Over a wide cutoff 
range the EFT prediction reproduces the LQCD result within its error, evidence that the 
EFT in LO captures the essence of the strong-interaction dynamics. As D varies (at hxed 
A) within the error bars of R 3 , B 4 also changes within its error bars. The estimate of a 
40% error in LO EFT is likely conservative, indeed. 

The power of EFT is the relative ease with which it can be solved for more-body 
systems. Binding energies for A > 5 are predictions that extend LQCD into new territory. 
Using A = 2fm“^, the authors of Ref. [I] searched unsuccessfully for excited states in 
A = 2,3,4 systems. Similarly, they found no evidence of ^n droplets, for which the 
ground-state binding energy coincided with the two-body threshold. Results [1] for the 
A = 5,6 ground states at A = 2 fm“^ are shown in Tab. [T], with errors estimated from 
the EFT truncation. For ®He a bound state with binding energy B^ = 98.2 MeV for 
A = 2 fm~^ coincided with the four-body threshold for A = 4fm~^. The ®Li ground state 
for A = 2 fm~^ was found at Rg ~ 122 MeV. In this case the error in K^ax extrapolation 
was about 3 MeV, which is somewhat larger than for lighter systems but still small 
compared with input and truncation errors. Calculations with AFDMC at larger A are 
in progress [ 6 |. 

These results are afflicted by considerable error bars. Even though the 40% assigned to 
the EFT expansion is likely an overestimate, the LQCD input itself has large uncertainties 
of about 25%. With this caveat, the trend of the results is surprising. There is a qualitative 
difference with A = 2 at the physical pion mass, because the dineutron is bound at larger 
masses. (For the effects of the dineutron scattering length on light nuclei in Pionless EFT 
at physical pion mass, see Ref. [S].) But this is consistent with other binding energies, 
which are all larger, and all larger by roughly similar amounts. The gap at A = 5, familiar 
in nature, seems to survive the increase in pion mass. And Rg/S ~ 20 MeV, similar to 
lattice ^He for which R 4/4 25 MeV. This suggests that nuclear saturation might not 

be tremendously sensitive to the pion mass. Overall, it seems that the lattice world at 
= 805 MeV is not that different from our own, with Ba/A scaled by a factor 4 or 5. 


5 Outlook 

Natural light (A < 4) nuclei are halo-type systems in the sense that they have sizes 
large compared the range of the force, and are thus described by Pionless EFT. One 
might expect this feature to result from hne-tuning, and to hnd very different worlds 
at unphysical pion masses. Surprisingly, the hrst LQCD calculations and their EFT 
extrapolations seem to suggest the opposite. Perhaps pions do not play as decisive a role 
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in low-energy QCD as we are used to think, and some of the dehning properties of nuclei 
are relatively insensitive to the value of the pion mass. If this is true, it should have 
implications for the use of nuclei in tests of the variability of fundamental constants |48] . 

Of course, at this point these are only hints. The same exercise can be, and is being |5], 
carried out in LO with the LQCD data at = 500 MeV [7j. EFT extrapolations can be 
repeated at other pion masses as LQCD results appear. More urgently, EFT calculations 
need to be performed for A > 4 at higher cutoffs to conhrm renormalizability, even at 
the physical pion mass. The scattering lengths and effective ranges at = 800 MeV 
[9] give just enough input for an NLO analysis, which should allow stronger statements 
about the convergence of Pionless EFT at unphysical pion masses. Finally, at all values 
of pion mass, one should increase A [6] to conhrm trends in Ba/A and to seek the limit 
of applicability of the EFT. 

On a longer time frame, as pion masses in LQCD drop sufficiently, one can use Chiral 
EFT to extrapolate further down in pion mass and fully predict real nuclei, following the 
same steps as in Refs. [lElE]. It is a lot to do, but it promises to fulhll a longstanding 
dream of nuclear effective held theorists, teaching us much about the connection between 
QCD and nuclear physics. 
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